The effect of finite container size on granular jet formation 
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When an object is dropped into a bed of fine, loosely packed sand, a surprisingly energetic jet 
shoots out of the bed. In this work we study the effect that boundaries have on the granular jet 
formation. We did this by (i) decreasing the depth of the sand bed and (ii) reducing the container 
diameter to only a few ball diameters. These confinements change the behavior of the ball inside 
the bed, the void collapse, and the resulting jet height and shape. We map the parameter space 
of impact with Froude number, ambient pressure, and container dimensions as parameters. From 
these results we propose a new explanation for the thick-thin structure of the jet reported by several 
groups [J.R.Royer et al, Nature Phys. 1, 164 (2005)], [G.Caballero et al, Phys. Rev. Lett. 99, 
018001 (2007)] and [J.O. Marston, et al, Physics of Fluids 20, 023301 (2008)]. 

PACS numbers: 45.70.-n, 47.55.Lm, 47.57.Gc 



I. INTRODUCTION 

Granular materials consist of discrete particles which 
interact mainly through contact forces. In large quanti- 
ties they can behave like a solid, a liquid, or a gas but 
often behave differently from what would be expected of 
these phases A marked example is the impact of an 
object into a bed of sand. When dry air is blown through 
such a bed all contact forces between the individual par- 
ticles are broken and after slowly turning off the air flow, 
the bed settles into an extremely loosely packed solid-like 
state. When a ball is dropped in such a bed, one observes 
a splash and a jet, strikingly similar to the ones that are 
seen when the same object is dropped into a liquid. 

Research interest in this granular jet started when 
S.T. Thoroddsen and A.Q. Shen first reported this phe- 
nomenon in 2001 0, in a study with the objective to gain 
insight into the importance of surface tension on jetting 
in general and the properties of flowing granular materi- 
als. Since these results, several aspects of the formation 
of the granular jet have been studied. The influence of 
the impact velocity onto the jet height for impacts on 
a bed of very loose sand was investigated in Q ■ Using 
a pseudo two-dimensional setup, numerical simulations 
and comparisons to water impact experiments, a model 
for the jet formation was proposed that is based on cav- 
ity collapse: The impacting ball creates a cavity in the 
sand bed which collapses due to the hydrostatic pressure 
in the sand and leads to two vertical jets. One jet is ob- 
servable above the bed and the other one is going down 
into the bed The series of events is concluded by 
a "granular eruption" at the surface of the sand which 
was attributed to the surfacing of an air bubble that is 
entrapped during the collapse. 

The influence of the ambient pressure on the forma- 
tion of a granular jet was first studied by Royer et al. [1]. 
They observed that at lower ambient pressures the jet 
reaches less high and also reported a puzzling thick-thin 



structure at lower pressures. Using X-ray radiographic 
measurements, they were able to look inside the bed and 
then proposed the following mechanism to explain this 
structure: the thick jet is caused by the compressed air 
in the cavity pushing up bed material, forming the thick 
part of the jet [1, ID, [||. The thin jet was attributed 
to the hydrostatic collapse as formulated in ||. Sub- 
sequently, the thick-thin structure was also observed at 
atmospheric pressure by increasing the ball size in the 
same container, which suggests -in contrast to the ear- 
lier explanation- that the structure may be a boundary 
effect [3] • Marston et al also found a thick-thin structure 
by decreasing the packing fraction, and they too found 
that this effect is more pronounced for a larger ball 
It is the exploration of the formation of this thick-thin 
structure that constitutes the main motivation for the 
work described in the current paper. 

In parallel to the research concerning the formation 
of the granular jet, quite some effort was made to un- 
derstand the motion of an object moving through a 
granular medium. Different drag force laws were pro- 
posed d, [l(| [0, H, [H, [3, [H[ ! culminating in a model 
containing a hydrostatic term that linearly depends on 
the depth inside the bed and a dynamic term which 
is proportional to the square of the velocity of the ob- 
ject [H, [H|. The influence of the ambient air pressure 
on this trajectory was investigated in 0] where it was 
shown that the drag force reduces at high ambient pres- 
sure. Another important issue is the interaction between 
the impacting ball and the container boundaries. Nelson 
et al. found that "the presence of sidewall causes less 
penetration and an effective repulsion" [l|| [TtJ • 

In this paper, we present experiments in which the 
size of the container has been systematically reduced. 
We did this by (i) decreasing the depth of the sand bed 
(section IIV[) and (ii) reducing the container diameter to 
only a few ball diameters (section [Vj . We explore how 
these confinements change the behavior of the ball inside 
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the bed, the void collapse, the resulting jet height and 
shape, and the presence of a granular eruption, which 
was only observed in part of the parameter space cov- 
ered in this study. All of the observed phenomena are 
explained within the context of a simple hydrostatic col- 
lapse model Q together with a drag law for the trajectory 
of the ball inside the sand [l5[. Finally, we propose an 
explanation for the presence of an eruption and a new 
mechanism for the thick-thin structure reported by sev- 
eral groups mentioned above. 

The paper is organized as follows: In Section[TT]wc start 
with the introduction of the drag law and the hydrostatic 
collapse model that lie at the heart of the analysis of 
this paper. Subsequently we discuss our experimental 
setup in Section ITTTl after which we present our results for 
impacts in confined settings. In Section IIVI we discuss 
the influence of the proximity of the bottom, after which 
we turn to the influence of the side walls in Section fVl 
Finally, in Section [VI] we discuss the thick-thin structure 
and end with conclusions in Section [VIH 



II. DRAG LAW AND HYDROSTATIC 
COLLAPSE MODEL 
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FIG. 1: Schematic representation of the impact of a ball into 
a sand bed, indicating the time and length scales that play 
an important role in the analysis of the experimental work in 
this paper, as described in the main text. 



In this Section we review the drag law and the 
Rayleigh-type collapse model that constitute the ingre- 
dients of the hydrostatic collapse model first introduced 
in Q and form the theoretical backbone against which 
the experiments will be analyzed. Before doing so let 
us recall the succession of events after an impact of a 
sphere in a granular bed, presented in Fig. [TJ which in- 
volves the introduction of several time and length scales 
that are crucial to the analysis in the following Sections. 
At a time t = 0, the sphere impacts on the granular 
medium with a velocity vq. A splash is created and the 
ball penetrates into the sand bed. The void created by 
the ball collapses in a time t c (closure time) and leads to 
the formation of two jets (one of them is visible above 



the bed and the other one is going down into the bed). 
The closure depth -also known as the pinch-off point- is 
denoted as z c and the position of the ball inside the sand 
at that time as z(t c ). At this time the ball will typically 
still be moving inside the sand bed, until, after a time t s , 
the ball reaches its final depth zf and stops. Finally, the 
entrapped air bubble slowly rises inside the sand bed and 
leads to a granular eruption at the surface at t = t einp . 

The first ingredient of the hydrostatic collapse model 
concerns the motion of the ball with diameter d through 
the sand bed. To describe the trajectory of the ball (z(t) 
is the depth of the ball at a time t), we use the law 
introduced by Tsimring [l4[ and Katsuragi [lj|. The 
drag force is decomposed into two terms: The first one, 
the hydrostatic term, involves Coulomb friction as well 
as the force needed to displace material against the hy- 
drostatic pressure and is proportional to the depth and 
was introduced in this context in 
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We here write 
-^hydrostatic = kz where k is a constant. The second 
term is an quadratic drag independent of the depth, 
^dynamic = otv 2 , interpreted as the quadratic force re- 
quired for the projectile to mobilize a volume of granular 
media with density p g proportional to the volume of the 
ball [22^ . Adding gravity, this leads to the equation of 
motion: 



mg 



a v 



(1) 



with initial conditions z(0) = and i(0) = vq. 

The second ingredient regards the dynamics of the hy- 
drostatic collapse of the void that is formed by the ball. 
The radius of the void at a time t and a depth z, R(z, t), 
evolves from the two-dimensional Rayleigh-type equa- 
tion, in which, for each depth z, the collapse is driven 
by the hydrostatic pressure p g gz at that depth 0] 



(RR + R 2 ) log 



R 



\r 2 = gz , 



(2) 



where R denotes the time derivative and i?oo is a constant 
of the order of the system size. Under the assumption 
that the cavity that is created is approximately cylindri- 
cal, i.e., with the same diameter (d) as the impacting ball, 
the initial conditions are R(0) = d/2 and R(0) = 0. By 
rescaling lengths with the ball radius d/2 and time with 

d/{2^J~gz) (i.e., R = 2R/d, R = Rj \f~gz, etc., where the 
dot on a dimensionless variable denotes a derivative with 
respect to dimensionless time) , Eq. ^ can be written in 
dimensionless form 



~~ ~ 2 R 
{RR + R ) log — 

Roo 



2 R 



(3) 



together with initial conditions R(Q) = 1 and R(0) = 
0. With these initial conditions this equation has a 
unique solution R(t), from which we obtain a constant 
dimensionless collapse time i co ii. It now follows im- 
mediately that the (dimensional) collapse time t co n [= 
tcottd/(2y/gz)] scales as ~ dj(2J~gz). 
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Finally, we can combine the above two ingredients to 
determine the position and the time of closure. The total 
time that elapses from the impact to the collapse of the 
cavity at any depth z is given by: 

*tot(z) = t pass {z) + tcou(z) ■ (4) 

where t paS s is the amount of time the ball takes to reach 
depth z (obtained from solving the drag law) and t co n 
is the time needed for the collapse at a depth z. The 
closure depth is the depth which minimizes equation |4] 
The closure time corresponds to the total time at the 
closure depth (t c = t to t{z c ))- 

III. EXPERIMENTAL SETUP 

In the previous section, we have introduced the the- 
oretical framework for the analysis of the phenomenon. 
We now turn to the description of the experimental setup, 
which is sketched in Fig.[2j It consists of a container with 
a height of 1 m and a square cross section of 14 x 14 cm 2 , 
which is filled with sand grains, nonspherical and slightly 
polydisperse in size (between 20 and 60 /jm); the den- 
sity of the granular medium is 2.21 g/cm 3 and its angle 
of repose 26°. As described in Q, before each experi- 
ment, the sand is fluidized by blowing pressurized dry 
air through a sintered plate at the container bottom. Af- 
ter slowly turning off the air flow, the bed reproducibly 
settles into a static, loose, weakened state (volume frac- 
tion 41 %). The airtight system can be slowly evacu- 
ated to perform experiments at lower ambient pressures 
p (the pump speed is low enough not to irreversibly al- 
ter the packing fraction). Then a steel ball of diameter 
d = 1.6 cm and mass m = 16.5 g is dropped into the 
sand from different heights H which controls the impact 
velocity vq = \J2gH, where g is the acceleration of grav- 
ity. Thus, the impactor is characterized by a single di- 
mensionless number, the Froude number (Fr), defined as 
Fr = 2vl/{gd) = AH/d. 

The impact is recorded with a high-speed camera 
(Photron Ultima APX-RS) at 1500 frames per second. 
For the measurements a uniform lighting from behind is 
needed to obtain better movies with higher contrast be- 
tween the objects and the background. This is realized by 
positioning two light sources and a diffusing plate behind 
the container. 

In order to obtain the trajectory of the sphere in- 
side the sand, we attach a wire with markers which re- 
main above the sand during impact and are imaged with 
the high-speed camera. This procedure is explained in 
greater detail in Section fVl-A. 

We use two ways to confine the impact and jet forma- 
tion experiment. First of all, we study the influence of the 
bottom of the container by reducing the height at which 
the container is filled with sand down to a few ball diam- 
eters (Scction llV|) . Second, to investigate the influence of 
the closeness of the side walls, we insert PVC cylinders of 
varying diameters into the sand, such that the cylinder 



axis coincides with the trajectory of the ball inside the 
sand. In this procedure sufficient care was taken that the 
presence of the cylinder did not disturb the fluidization 
and settling process of the sand bed (Section [Vj. 
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FIG. 2: Setup: (a) perspex container, 14 x 14 x 100 cm 3 , (b) 
pneumatic release mechanism, (c) Photron Ultima APX-RS, 
(d) two light sources with diffusing plate, (e) pressurized, dry 
air source, (f) computer and (g) vacuum pump with pressure 
gauges. 



Time and position of collapse 

When traveling through the sand bed, the ball creates a 
cavity. The shape of the cavity is obtained using a pro- 
filometer similar to the one described in [l8[ (see Fig. [3]). 
A diode laser sheet with wavelength of 680 nm strikes the 
granular media at an angle 9, fixed arbitrarily at 55°. Us- 
ing a mirror and a high-speed camera, we can measure 
the horizontal projection of the points where the laser 
sheet touches the sand from above. When the surface is 
flat, this projection is a straight line parallel to the y- 
direction; the coordinate of a point on this straight line 
is (xi, y). When the surface is perturbed, the projection 
appears to be a curved line. For any point on this line 
with coordinate (x(y), y)) the depth of the surface can 
be deducted -as a function of y- from xi and x{y) 

z(y) = (x(y)-x l )t;m(6). (5) 

If we assume rotational symmetry of the cavity around 
the center of the ball [denoted as (x c ,y c )] we can in ad- 
dition deduce the radius of the cavity at all these depths 



R(z{y)) = ^(x{y)-x c y + {y-y c y. (6) 

By analyzing each of the high speed imaging recordings 
in this way, we can obtain the cavity profile R(z,t) as 
a function of both depth z and time t (up to a certain 
maximum depth that is set by the laser sheet angle 9). 
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FIG. 3: (I) Laser profile-meter. A diode laser sheet (a) is 
directed onto the surface at an angle 9. Using a mirror (b) and 
a high-speed camera (c), images of the surface are recorded. 
(II) Schematic view of the resulting surface. The dashed line 
represents the laser sheet when the surface is flat and the 
continuous line the laser sheet when the surface is perturbed. 
The local deviation Sx — x(y) — xi of the laser sheet is related 
to the vertical coordinate Sz = z(y) of the surface, (d) is the 
center of the cavity, from which the cavity radius R(z) can be 
deduced. 



IV. INFLUENCE OF THE BOTTOM: A 
SHALLOW BED 

Now that we have introduced the experimental setup, 
we will continue with the discussion of our results: In this 
Section we start with what is observed in a less-filled con- 
tainer (i.e., a shallow sand bed) and in the next Section 
proceed with the discussion of what happens when the 
diameter of the system is decreased. 

Before turning to the case in which the proximity of the 
container bottom becomes important, let us first recall in 
table [J the results obtained in the usual unconfined case, 
here at Fr = 70 and ambient pressure: the container is 
large enough (D = 14 cm S> d) to avoid any effect of 
the surrounding walls and deep enough (the height of 
the sand bed, hbed is around 30 cm, that is 18.75 d) such 
that the bottom has no influence. 



Final depth zi 


Stopping time t Bt o P 


collapse time t c 


lid 


110 ms 


51 ms 


closure depth z c 


jet height hj e t 


eruption time t e rup 


2d 


18.5 d 


510 ms 



TABLE I: Results obtained at Fr = 70 and p = 1 bar in the 
usual unconfined case, i.e., in a deep bed with ftbed = 18.8 d. 
These values will be used as reference values in the discussion 
of the experimental results. 

We modified the height of the sand bed, hbed by simply 
adjusting the amount of sand in the container. The first 
and most conspicuous effect is that below a certain depth 
of the sand bed the impacting sphere is stopped abruptly 
by its collision with the container bottom, rather than 
slowly being stopped by drag as happens in the uncon- 



fined case. In this way, decreasing the depth of the sand 
bed allows us to look at the influence of the final depth 
of the ball, Zf, and the cavity size on the jet and the 
eruption. 



A. Influence on the jet 




h b ed / d = 5 h be[) / d = 4.1 h bed / d = 2.6 h bed / d = 2.1 



FIG. 4: Influence of the height of the sand bed hbed on the 
shape and height of the jet for Fr = 70 and p = 1 bar: Images 
of the jet, taken at 0.12 s after the ball impact for four different 
bed heights, decreasing from left to right. Below a threshold 
there is a clear change in height and width of the jet. 

In Fig. 31 we show four images from the jet formed 
when the ball is dropped into the sand bed for Fr = 70 
and ambient pressure. While reducing hbed, there is 
no change in the jet shape or height down to a certain 
threshold. Below this threshold, the jet reaches less high 
and becomes broader, most notably at the top. The max- 
imum height of the jet, h ie t, is obtained by measuring the 
initial jet velocity i>j et as soon as it appears above the sur- 
face of the sand (using energetic arguments, h$ e t cx v? et ). 
The initial jet velocity Vj e t is plotted as a function of hbed 
in Fig. [5j For hbed higher than 3 d, v- }e t saturates to its 
undisturbed value of approximately 3.2 m/s. Reducing 
hbed below 3 d, v- 3e t decreases rapidly. When we reduce 
the ambient pressure to p — 100 mbar, we find the same 
behavior (see Fig. [5j) although the crossover takes place at 
a slightly higher value of hbed- Remarkably, in both cases 
this decrease does not happen at the depth at which the 
ball is stopped by the bottom (which would be around 
^bed = 11 d and hbed = 6d for p = 1000 mbar and 
p = 100 mbar respectively) but at a much lower depth 
of hbed ~ 3d. 

This can be explained as follows: The closure depth, z c , 
remains unaltered by the presence of the bottom (which 
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FIG. 5: Initial velocity of the jet, Vj at as a function of the 
height of the sand bed, hbed for Fr ~ 70 and p = 1000 mbar 
(+) and p = 100 mbar (o). There is a sharp threshold below 
which the initial jet velocity rapidly decreases. The dashed 
lines represent the undisturbed values of Vj et , measured in a 
deep bed (/ibcd = 18.8 d). 



below /ibcd = Hd only makes the ball stop earlier and 
less deep) until the bed depth becomes less than the sum 
between the position of the unconfined collapse (2 d, see 
tabic fl} and the diameter of the ball. Below this value, 
the collapse happens on top of the ball leading to a less 
directional top of the jet which has a more or less spheri- 
cal shape; moreover the closure depth decreases when the 
bed becomes smaller and so does the initial jet velocity. 
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FIG. 6: The time t elup when the granular eruption at the sur- 
face starts is plotted as a function of the height of the sand 
bed hbed, for Fr = 70 and p — 1 bar (black open circles). 
Measurement points with t CIup — correspond to those cases 
where no eruption was observed. The experimental regimes 
with and without eruption are separated by the vertical black 
lines. The grey region represents the region where no air bub- 
ble is entrapped. The thin blue and red lines represent the 
different time scales that are involved in the problem t\ is the 
time the air bubble needs to reach the surface (black thin line) 
and t2 is the time the air bubble needs to diffuse within the 
sand bed (red thin line). When t\ is smaller than £2, an erup- 
tion is expected; this is depicted by the continuous thick blue 
line. The different regions obtained from the timescale ar- 
gument qualitatively correspond to the experimental results. 
More details about the way in which t\ and t2 are estimated 
are provided in the main text and in footnote [2^ ]. 



B. Influence on the eruption 

Providing that the void collapse does not happen di- 
rectly at the ball, an air bubble is entrapped. The volume 
of this bubble can be estimated as: 

Vbubble oc h 2 ad (z(t c ) - z c ) ~ d 2 (z(t c ) - z c ) , (7) 

where z(t c ) is the position of the ball at closure and h ra d 
is the radial length scale of the bubble, which can be 
approximated by the diameter of the ball. The bubble 
slowly rises through the sand and can lead to a violent 
granular eruption. However, this eruption is not always 
observed. To study when and why this is the case, in 
Fig. [5] we plot the time between impact and eruption, 
iorup, as function of the height of the sand bed, hbed- 
(Note that measurement points with t clup = corre- 
spond to those cases where no eruption was observed.) 
Up to a certain threshold, which is around 4.8 d, no erup- 
tions can be observed. This can be attributed to the fact 
that, while rising, small air bubbles just dissolve into the 
sand bed before reaching the surface. When the bed gets 
deeper, the air bubble reaches a certain critical volume 
V*, above which a granular eruption can be seen. From 
the experimental results, this size found to be around 
V* ~ d 2 (z(t c ) — z c ) ~ 3.8d 3 . Then remarkably, above 9 d 
the eruptions disappear again and reappears only when 
hbed > 14 d. 



This peculiar behavior can be understood, at least 
qualitatively, from the competition of the two time scales 
corresponding to the two mechanisms the air in the bub- 
ble has to escape from the bed: 

• The bubble needs a time t\ to reach the surface. 
First of all, for hbed < 3 d, the collapse happens 
on top of the ball, and no air bubble is entrapped. 
Between 3d and 5.5 d, the position of the ball at 
closure, z(t c ), increases and so does the volume of 
the air bubble; in this region, t\ decreases. While 
increasing the sand depth even further, the volume 
of the air bubble remains constant, but the initial 
position of the bubble is deeper since the entrapped 
bubble will follow the ball until it stops. Therefore 
t\ will increase until hbed is equal to lid which is 
the final depth of the ball in the unconfined case. 
Above this value, there is no change on the final 
depth and t\ becomes independent of the depth of 
the sand bed. This is depicted by the thin black 
line in Fig. [5] More details about the way that t\ 
is estimated are given in [23| . 

• The air in the bubble trapped by the collapse es- 
capes -in the dissolution time t-r- through pressure 
driven flow through the porous bed. Factors that 
affect this process are the size of the bubble (which 
determines the amount of air that needs to escape), 
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the pressure of the air (which approximately equals 
the hydrostatic pressure in the surrounding sand), 
and the length of the path the air needs to travel. 
For this last quantity we need to consider that air 
can both escape through the top and through the 
bottom of the bed, the latter due to the presence 
of the sintered plate. These factors combine into 
Darcy's law: Q oc AP/H, where the flow rate Q 
depends on the pressure difference AP driving the 
flow and the path length H. Turning to Fig. [5] 
again, for /ib c d < 3 d, no air bubble is entrapped. 
Between 3 d and 5.5 d, z(t c ) -and therefore the bub- 
ble size- increases, leading to a steep increase of 
the dissolution time £2- Upon increasing the sand 
depth beyond 5.5 c? the bubble size remains con- 
stant but the pressure inside the bubble increases. 
From Darcy's law we thus find that decreases. 
This continues until /ibed is equal to lid beyond 
which the ball does not reach the bottom of the 
plate anymore. Note that until this point the short- 
est (and therefore chosen) path for the bubble to 
dissolve is towards the bottom of the container. If 
we now increase /ibed beyond 11c? this shortest path 
starts to grow, and with the path, using Darcy's 
law, also the dissolution time. This is captured by 
the thin red line in Fig. [HI More details about the 
estimation of ti are provided in 23]. 

As a result, an eruption can be only seen if the time £2 
becomes larger than t\. This is expressed by the contin- 
uous thick blue line, in qualitative agreement with the 
experimental behavior. 



V. INFLUENCE OF THE SIDE WALLS 

In the previous section, we discussed the influence of 
the bottom of the cavity on the process of object penetra- 
tion and jet formation and found that, if the sand depth 
is fixed at 30 cm, there is no effect of the bottom on the 
jet formation process. Fixing this bed depth, we now 
turn to study the effects of the side walls of the container 
on the complete scries of events leading to the jet. For 
this, some cylinders of different diameters D are placed 
inside the sand during the fluidization process: we choose 
D = 4.2 cm, 6 cm, 8.5 cm, 10 cm and 12.5 cm. In this 
way we change the aspect ratio D/d from 2.6 to 7.8. 



A. Ball trajectory 

The first thing that happens upon impact of the ball 
onto the surface is that it penetrates and creates a void 
inside the sand bed. The question we ask in the next sub- 
section concerns the influence of the container diameter 
on the drag force experienced by the ball during its mo- 
tion through the bed. To monitor the trajectory of the 
ball, a wire with a tracer is attached to the ball. Using 



a high-speed camera (1500 frames per second) and im- 
age analysis, we obtain the trajectory of the tracer and 
therefore the trajectory of the ball z{t). To keep the wire 
tense an extra friction device and a light counterweight 
are used , which have the effect that the ball experiences 
a downward acceleration due to gravity which is approx- 
imately 10 % smaller than g. The actual acceleration 
is measured during the "free fall" part of the trajectory, 
and the results presented here have been corrected for 
this effect. 

In the top two plots of Fig. we compare the trajec- 
tories of the ball at ambient pressure for an impact with 
Fr = 25 and for two diameters of the confining cylinder 
(D = 6.0 and 10.0 cm). We can fit the experimental tra- 
jectories using the model introduced in Section HT1 fEq . [Tj) 
using a and k as fitting parameters. The agreement be- 
tween the model and the experiments is very good (see 

Fig. uj) . 

Decreasing the diameter of the container surprisingly 
increases both the final depth of the ball, Zf and the time 
to reach the final depth, t s . In Fig. [5^, we report the 
final depth Zf as a function of the container diameter 
at different pressures for Fr = 25. There is a clear de- 
pendance: The final position of the ball is deeper for a 
smaller container. Also, the influence of the boundaries 
for this Froude number is less pronounced at small pres- 
sures. We conclude that for Fr = 25 the drag force the 
ball experiences becomes smaller for small containers. 

But what happens at higher Froude numbers? In 
Fig. [5b, we report the final depth, Zf as a function of the 
container diameter for Fr = 75. At first glance the behav- 
ior now seems completely opposite to what we observe at 
small Froude number, as the final depth now decreases 
with decreasing container diameter: To be more precise, 
at atmospheric pressure the final depth stays more or 
less constant and at lower pressures there is a decrease in 
Zf with decreasing container diameter. So now the drag 
force seems to be larger for small container diameters. 

To understand this difference, we have to separately 
look at the behavior of the hydrostatic and quadratic 
drag force: After all, for small Froude numbers we ex- 
pect that the hydrostatic drag —kz will dominate the 
behavior of the ball, whereas for higher impact velocities 
it is expected that the quadratic drag av 2 will start to be- 
come increasingly more important during the motion of 
the ball. To this end, in Fig. [5] we plot k and a as a func- 
tion of container diameter for three different pressures. 
Each value represents the average parameters obtained 
from fits to the trajectory data analogous to the ones of 
Fig. 0over a range of Froude numbers from 25 to 100 [24(. 
As shown by Caballero 0] , the hydrostatic force depends 
on the ambient air pressure: k decreases with p roughly 
as p -1 / 2 . Our findings are consistent with this observa- 
tion, also for smaller container diameters (not shown). 
Next to this we find that k increases quite steeply with 
D, which is consistent with the lower drag experienced 
by the impacting ball at smaller container diameters at 
low Froude numbers. Physically, the behavior of the hy- 
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FIG. 7: (top) Depth of the ball z(t) (a) and its velocity v(t) (b) as a function of time after impact for Fr = 25, p = 1000 mbar, 
D = 10 cm (o) and D = 6 cm (□). The lines correspond to a fit using equation Q] with k = 4.525 N/m and a = 0.132 kg/m for 
D = 10 cm and k — 1.695 N/m and a = 0.118 kg/m for D = 6 cm. (bottom) Depth of the ball z(i) (c) and its velocity v(t) (d) 
as a function of time after impact t for Fr = 75, p = 50 mbar D = 4.2 cm (o) and D = 12.5 cm (□). Again, the lines correspond 
to a fit using equation [1] with k — 14 N/m and a = 0.281 kg/m for D = 4.2 cm and with k = 13.5 N/m and a = 0.111 kg/m 
for D = 12.5 cm. Within the smallest container, and only at low pressure, we observe anomalous behavior: The ball reaches a 
plateau in which the velocity remains constant before going to zero again at larger times. Clearly, the model fails to describe 
the data in this case. 



drostatic drag force can be understood using a similar 
argument as p} : When the ball passes through the sand, 
an air flow is created around it which fluidizes the sand 
bed and reduces the drag force. This effect is expected 
not only to be more important at higher pressure but 
also when the container diameter becomes smaller: Near 
the wall, the velocity of the interstitial air is required to 
be zero and, since the same amount of air needs to be 
displaced, the flow will be more important if the aspect 
ratio D/d is small. Consequently, the hydrostatic drag 
force will be lower. 

Figure [5b contains the coefficient a of the quadratic 
drag term av 2 . Clearly, a becomes larger for smaller 
container diameters but the difference is hardly as pro- 
nounced as was the case for k. This accounts for the ob- 
servation that at some point, for larger Froude number, 
the drag does become larger when the container diameter 
is decreased. 

Finally, in Figs. [5^ and b there is one exceptional value: 
For the smallest container diameter (D/d = 2.6) and the 
lowest pressure (50 mbar) the fitted values of k and a 
turn out to strongly depend on the Froude number. The 
bottom two plots in Fig. which contain two trajec- 
tories at 50 mbar for the largest and the smallest con- 
tainer diameter, reveal the reason why: Whereas for the 
biggest container (D = 12.5 cm), the behavior is similar 
to the behavior described for Fr = 25, for the smallest 



one (D = 4.2 cm) it is qualitatively different. Whereas 
the agreement between the experimental and the com- 
puted trajectory sill seems to be reasonable (Fig. [7J;), 
the velocity of the ball (Fig. [7ji) presents large discrep- 
ancies: The measured ball velocity doesn't decrease to 
zero gradually, but first slows down until it reaches a 
plateau at constant velocity where it stays for a while 
before slowing down until it stops. This behavior is iden- 
tical to the one observed in the X-ray experiments of 
Royer et al. 0] , in which the container needed to be kept 
small. That behavior therefore is likely to be a boundary 
effect. We believe that the origin of the plateau lies in 
a depth-independent force between the ball and the wall 
(mediated by force chains) which is dominant over the 
hydrostatic drag force and, together with the quadratic 
drag force, balances gravity at the plateau velocity. At 
some depth, the Coulomb drag force takes over, slow- 
ing the ball down to zero. Obviously, the model cannot 
be valid in this situation and an extra force due to the 
ball/ wall interaction should be taken into account. 

B. Collapse of the cavity 

The second issue that we want to address in this Sec- 
tion is the influence of the container diameter on the 
collapse of the cavity. We study the dynamics of the col- 
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FIG. 8: Final depth Zf as a function of the container diameter 
D, at different pressures, for a) Fr = 25 and b) Fr = 75. The 
final depth is divided by the final depth for the unconfmcd 
case in order to emphasize the deviations due to the proximity 
of the boundaries. The dashed lines are a guide to the eye to 
separate the different pressures. 



lapse of the cavity at closure depth using the profilomet- 
ric method described in detail in Section IIIII In Fig. [TUl 
the radius of the cavity is plotted as a function of time 
t — ipass for two different diameters at atmospheric pres- 
sure where t pass is the time needed for the ball to reach 
the closure depth z c . We can clearly distinguish a slight 
expansion of the cavity followed by a strong contraction. 
The collapse accelerates towards the pinch-off. Due to 
experimental limitations we do not have enough spatial 
resolution to obtain data points close to the pinch-off. 
The void dynamics is in qualitative agreement with the 
behavior predicted by the 2D Rayleigh-Plesset equation 
described in section [IT] (Eq. O . Whereas the expansion 
turns out to be weak and more or less independent of 
the container diameter, the contraction and the closure 
strongly depend on it. A plausible explanation would be 
that, for small containers, less sand is involved in the col- 
lapse. Therefore, the hydrostatic pressure which drives 
the collapse is not as sustained as for a larger container, 
explaining why the collapse takes longer for a smaller 
container (Fig. [TU|) . 

In Fig. [11] we plot the closure depth z c and the closure 
time t c . We find that t c increases and z c decreases when 
decreasing the container diameter. This decrease of the 
closure depth is generic: Also for small Fr, where Zf ac- 
tually increases, we find a decrease of z c . The fact that 
a decrease of the closure depth z c implies an increase of 
the collapse time t co \\ can be understood from a reduc- 
tion of the driving pressure (oc gz c ) and the availability 
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FIG. 9: a) n and b) a as a function of cylinder diameter 
D for different pressures p. For almost all values of D and 
p variations of both k and a are within the measurement 
error and each point is obtained from an average over a range 
of Froude numbers from 25 to 100. Only for the smallest 
container (D/d — 2.6) and the lowest pressure (50 mbar), 
there is a strong dependance of k and a on the Froude number; 
the model is not valid in this situation. Plot b) reveals that 
for large Fr the quadratic drag takes over for small cylinder 
diameters leading to less intrusion of the ball (Fig. [Sb). 




FIG. 10: Dynamics of the cavity collapse at closure depth for 
two container diameters D — 4.2 cm (□) and D — 10 cm (o). 
Here, Fr = 70 and p = 1 bar. The time has been rescaled 
by multiplying with a factor 2^/gz c /d in order to show the 
results in a single plot. The continuous line correspond to a 
fit using the 2D Rayleigh-Plesset equation (Eq. 0. 
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above). 
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FIG. 11: (a) Closure depth z c as a function of the container 
diameter D for different pressures, (b) Closure time t c as a 
function of the container diameter D for different pressures. 
For all measurements Fr = 70. 
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FIG. 12: The jet height, hj e t as a function of the container 
diameter D for Fr = 25 (a) and Fr = 50 (b) at different am- 
bient pressures. The jet height is divided by the jet height 
in the unconfined case in order to see the deviations due to 
the proximity of the boundaries. For all pressures and Froudc 
numbers the jet height increases with increasing container di- 
ameter. The dashed lines are a guide to the eye to separate 
the measurement series at different pressures. 



C. Jet Height 



Now that we studied how the process of the formation 
and subsequent collapse of the cavity are influenced by 
the proximity of the side walls, we continue with the in- 
fluence of the diameter of the container D on the jet and, 
in particular, on the maximum height reached by the jet. 
In Fig.|T2l the jet height h ]ct is plotted as a function of the 
diameter for two Froude numbers (Fr = 25 and Fr = 50) 
at different values of the ambientpressure. Since it was 
already discussed extensively in [7| that the jet becomes 
smaller at reduced ambient pressure, we now divide hj e t 
by the jet height in the unconfined situation. We ob- 
serve that, while decreasing the container diameter, the 
jet reaches less high. This behavior is the combined re- 
sult of the reduction of the closure depth and the increase 
of the closure time with decreasing container diameter as 
was described in the previous subsection: The reduction 
of z c reduces the hydrostatic pressure that drives the col- 
lapse and the increase of the closure time is connected to 
the fact that -because of the reduced container diameter- 
there is less sand available during the collapse, making 
the driving pressure less sustained. Both factors con- 
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FIG. 13: The jet height hj ct as a function of the container di- 
ameter D for Froude = 100 and different ambient pressures. 
Again, there is a clear change in jet height as function of con- 
tainer diameter. Measurements at the highest Froude num- 
bers are not possible due to the surface seal (see text). The 
dashed lines are a guide to the eye to separate the measure- 
ment series at different pressures. 
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tribute to a decrease of the jet height. The rcscaling by 
the unconfined jet height also reveals that the influence 
of the boundaries is similar for all pressures and even for 
these two different Froude numbers. The unconfined be- 
havior is obtained when the diameter of the container is 
seven times higher than the diameter of the ball. 

At high Froude number (Fr = 100), the results can 
only be obtained for small containers, because, when the 
diameter is large, the jet collides with the splash which 
is being sucked into the cavity behind the ball at high 
ambient pressures. This is similar to the surface seal that 
has been observed for impacts on a liquid [2(| Hl[ . For 
this high Froude number the results are less conclusive, 
as can be seen in Fig. [T2J This is possibly due to the 
increased importance of the air flow caused by the ball 
when it is restricted to a smaller container diameter at 
these high impact velocities. 



D. Granular eruption 

Finally, we turn to the granular eruption that termi- 
nates the series of events. Since the container diameter 
has an influence on both the final depth and the closure 
depth, it is expected that the granular eruption will de- 
pend on the container diameter D. In Fig. 114a we report, 
for Fr = 100, the phase diagram indicating the presence 
of an eruption in (zf , D)-space. Note that the different 
values for the final depth Zf at fixed container diameter 
D have been obtained by varying the ambient pressure 
p. There is a marked dependancc on the container diam- 
eter D: More eruptions are observed in a small container 
than in a large container. 

This behavior can be explained using the influence of 
the side walls on the trajectory of the ball and on the 
collapse time together with the closure depth: For the 
same pressure, the closure time is larger, which leads 
to a deeper position of the ball at closure z(t c ), and at 
the same time the closure depth is smaller, increasing 
the size of the entrapped air bubble for small container 
diameters. If we replace the final depth in Fig. [T4k by 
the quantity (z(t c ) — z c )/d which is proportional to the 
volume of the entrapped air bubble [remember that it was 
argued that T4ubbic ~ d 2 (z(t c ) — z c ), see Eq. $7), which is 
subsequently divided by d 3 ], we remove the dependancc 
on D: In Fig. RBb the phase diagram is separated into 
two parts using a horizontal line representing a critical 
volume V* ~ 3.8c? 3 . This means that, independently of 
the diameter of the container, the bubble volume upon its 
formation has to be big enough to lead to an eruption. As 
was explained in Section HVl this is because the air bubble 
must have sufficient time to reach the surface before it 
has completely dissolved into the sand bed. Incidentally, 
the value for the critical volume determined from the 
phase diagram corresponds well to the value found in 
Section HVl 
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FIG. 14: (a) Phase diagram for the granular eruption at 
Fr = 100 as a function of the final depth Zf and the container 
diameter D. In both plots red open circles indicate parameter 
values where an eruption was absent, whereas blue plus signs 
stand for parameter values with an eruption. Note that Zf 
(which is a measured quantity) has been varied by using dif- 
ferent values of the ambient pressure p. (b) The same phase 
diagram, now as a function of the volume of the entrapped air 
bubble ((z(t c ) — z c )/d) and the container diameter D. The 
latter plot clearly indicates that the presence of the eruption 
is a function of the entrapped air bubble size only [25l ]. 



VI. JET SHAPE AND THICK-THIN 
STRUCTURE 

The proximity of the side walls and the bottom does 
not only affect the height of the jet but also its shape. 
One of the most prominent features is the thick-thin 
structure first described by Royer et al. 0, Q who studied 
the dependence of this structure on ambient pressure and 
Froude number. In the same work Royer et al. proposed 
a formation mechanism for the thick part of the jet based 
on the pressurized air bubble pushing sand into the thin 
jet originating from the pinch-off at the closure depth. 

In this Section we report, in addition to the Froude 
and pressure dependence, a pronounced dependence of 
the thick-thin structure on the proximity of the container 
boundaries. We propose an alternative model for the 
formation of the structure which semi-quantitatively ac- 
counts for the observed behavior of the phenomenon for 
the entire parameter space. 

In our experiments we can distinguish three different 
jet shapes, two of which exhibit a thick-thin structure: 
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(2) 




FIG. 15: Typical snapshots of the three distinct jet shapes 
observed in experiment: 1) Normal jet (for D = 10 cm, Fr = 
100 and p = 1000 mbar); 2) Thick-thin structure with sharp 
shoulder (for D = 8.5 cm, Fr = 100 and p = 100 mbar); 3) 
Thick-thin structure with a transition (for D — 10 cm, Fr = 
50 and p = 50 mbar). All snapshots show the fully developed 
shape of the jet at its maximum height. The snapshots are 
not on the same scale. 



(1) a 'normal' jet, in which the width of the jet gradually 
decreases from bottom to top, 

(2) a thick-thin structure with a sharp shoulder, where 
the thick lower part abruptly changes into a thin up- 
per part, 

(3) a thick-thin structure with a transition, characterized 
by a transient region in which the thick lower part 
gradually passes into the thin upper part. 

An example of each of the three jet shapes is shown in 
Fig. US 

First, we briefly look at the influence of the bed depth 
on these structures for a moderate Froude number (Fr = 
70). At atmospheric pressure we observe a 'normal' jet 
for all values of the bed depth /ibed (Fig. To observe 
a thick-thin structure we need to go to smaller ambi- 
ent pressures: At 100 mbar, a thick-thin structure with 
sharp shoulders can be observed in the unconfincd case, 
i.e., for large fibed (Fig. [l"5"|). Below a certain threshold 
(^bed < 4d), the thick-thin structure gradually disap- 
pears (Fig. IT5|) . This disappearance coincides with the 
disappearance of the entrapped air bubble below 3 d in 
which case the collapse happens more or less on top of 
the ball. 

The effect of the proximity of the side walls (within a 
sufficiently deep bed) is reported in the three phase di- 
agrams of Fig. [T71 where the jet shapes are classified as 



FIG. 16: Three snapshots of the shape of the jet at different 
values of the height htcd of the sand bed, taken 120 ms after 
the ball impact for P = 100 mbar and Fr = 70. For ftbed = 
4.1 d there is a clear thick-thin structure (with a transition 
region), which gradually disappears when the bed height is 
decreased to 3.4 d and 2.6 d. 



a function of container diameter and final depth [26| , for 
three different Froude numbers. For the lowest Froude 
number (Fr = 25), a thick-thin structure with a transi- 
tion is found only for the smallest Zf (which corresponds 
to the lowest pressure, p = 50 mbar) at intermediate 
container diameter. Thick-thin structures with a sharp 
shoulder are not found for this Froude number. When 
we increase the Froude number, the thick-thin-structure 
region is found to grow. Within the region contain- 
ing the transition variety of the thick-thin structure we 
observe the formation and growth of a region contain- 
ing the sharp-shoulder variety. Although the thick-thin- 
structure region grows to include the largest container di- 
ameters that we have used in our experiment [27| , thick- 
thin structures are never found in the smallest container 
diameter for the parameter space explored in this study. 

Now, which mechanism causes these structures? To 
answer this question we hypothesize a second collapse 
that happens on top of the ball forming a second jet. 
Such a second collapse can be motivated from experi- 
ments in a quasi-twodimensional setup [il^ and from X- 
ray measurements [3, H| , where multiple collapses have 
also been observed. The idea is as follows: Since the sec- 
ond collapse happens at a later point in time, the first jet 
is already well on its way in the formation process when 
the second one is being formed. Now, if the second jet 
can catch up with the first fast enough, it will hit its base 
and produce a thick-thin structure. When the time span 
between the two jets is too long however, the first jet will 
have (almost) fully formed and the collision of the second 
jet with its base will not disturb its shape [28|. 

To quantify this idea we need to estimate the interval 
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FIG. 17: Phase diagram of the observed jet shapes as a func- 
tion of container diameter D and final depth Zf for three dif- 
ferent Froude numbers: (a) Fr = 25; (b) Fr = 50; and (c) 
Fr = 100. The different final depths Zf (at fixed D) corre- 
spond to different ambient pressures (p = 50, 100, 200, 500, 
and 1000 mbar). The dashed lines are a guide to the eye to 
separate the different regions in the phase diagrams. 
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FIG. 18: Schematic drawing of the mechanism leading to the 
thick-thin structure. In case (a), the second collapse happens 
before a certain threshold time, such that the thickness of 
the layer of sand from the first collapse still is thin enough 
to be pushed up by the second jet and a thick-thin structure 
emerges. In case (b) we are above the threshold: The second 
jet collides with a thick layer of sand and is unable to disturb 
the formation of the first jet. 



between the time that the first jet is formed at the closure 
depth z c and the moment that the second jet reaches z c . 
This interval consists of the difference between the two 
closure times (t c ,2 ~ t c ) (where t C) 2 is the closure time of 
the lower collapse), summed with the time the second jet 
needs to reach z c , i.e., (zf — z c )/v2 with vi the velocity of 
the second jet. If this time interval is shorter than some 
threshold value T, we obtain a thick-thin structure, as 
visualized in Fig. [181 This leads to: 
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Let us illustrate the workings of this mechanism in an 
example: For Fr = 75 and p = 50 mbar we start from 
the largest container size where a thick-thin structure is 
visible. When decreasing the size of the container, the 
closure depth z c and the final depth Zf decrease following 
approximately the same behavior, such that the distance 
between the two collapses is more or less constant. Be- 
cause z f decreases, the hydrostatic pressure and therefore 
the velocity of the second jet decrease as well, such that 
the first term in Eq. JS]) increases. The same holds for the 
second term, because the closure time is found to increase 
with decreasing container diameter (cf. Fig lllb ). Thus, 
the left hand side of Eq. (J5J) increase with decreasing the 
container diameter, explaining why below a certain di- 
ameter the thick-thin structure disappears. 

To check wether the argument of Eq. j8]) also works 
quantitatively we approximate the several terms in the 
equation with experimentally known quantities. In the 
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first term V2 is proportional to the square root of the driv- 
ing hydrostatic pressure at depth Zf, i.e., v 2 = C^fgzf, 
with C constant. Because, similarly, for the velocity 
of the first jet we have Uj e t = Cy 'gz c , we find vi ks 
\fzifz c Ujet which is inserted into the first term of Eq. (JSJ). 
In turn, uj e t can be deduced from the jet height h- }at as 

In the second term, the unknown quantity is the sec- 
ond closure time t c ^2 -i.e., of the cavity just above the 
ball- which consists of the sum of the time t s the ball 
needs to come to a standstill and the time i co ii.2 the 
cavity needs to collapse at that point. Since accord- 
ing to the Rayleigh model discussed in Section |TT] the 
collapse times should scale as £ C oii.2 = C'd/ (2^/gzf) and 
icoli = C'd/ (2^/gz c ) respectively (with C constant), we 
have icoii.2 ~ \/Zc/ z f t co i\. Inserting all of the above in 
Eq. iJHJ) we obtain 
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thin structure. This is in agreement with the above 
mechanism: If an eruption is observed, this means that 
a relatively large air bubble must have been entrapped. 
This concurs with a large distance between the first and 
the second collapse point, which makes it unlikely that 
a thick-thin structure will be formed. Conversely, if a 
thick-thin structure is observed, this means that a (rel- 
atively small) air bubble must have been pierced by the 
second jet, which will facilitate its dissolution in the sand. 

Incidentally, our observations dismiss the mechanism 
for the formation of the thick-thin structure put forward 
in Q, where it was claimed that the rising air bubble 
pushes up bed material that forms the thick part of the 
jet. Contrary to the observations, such a mechanism 
would be stronger for a larger air bubble and cannot ex- 
plain why a thick-thin structure and a granular eruption 
cannot be seen at the same time. Next to this -for vary- 
ing container diameter- we observe both thick-thin struc- 
tures and normal jets for the same amount of entrapped 
air. 

Finally, based upon the present experiments we cannot 
exclude that the second collapse may need some down- 
ward motion of the sand bed that would be reinforced 
in a confined bed, but may become too weak to produce 
a jet in an unconfined bed. This would mean that if 
it were possible to increase the container diameter even 
further, the thick-thin structure may disappear again, as 
the phase diagrams in Fig.[l7|seem to suggest. This issue 
needs to be settled in future research. 



VII. CONCLUSIONS 
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FIG. 19: Phase diagram with on the vertical axis the left 
hand side of Eq. © and on the horizontal axis the container 
diameter D. The plot contains all measurements from Fig. 1171 
Short, black dashes indicate normal jets, intermediate, blue 
dashes the thick-thin structure with a transition, and long, 
red dashes thick-thin structures with a shoulder. 

In Fig. [19] we find a phase diagram in which all mea- 
surements from Fig. [T7| are plotted again, but now with 
the left hand side of Eq. © on the vertical axis. Clearly, 
all thick-thin structures (intermediate and large dashes) 
lie below some time-threshold, in agreement with the for- 
mation mechanism discussed above. The smallest con- 
tainer diameter forms an exception, in the sense that 
here thick-thin structures are also not found for time 
scales where they could have been expected (i.e., that 
lie clearly below the threshold T). This behavior maybe 
due to the fact that lack of material to sustain the col- 
lapse leads to an underestimation of the actual times in 
Eq. ([9]) . But in general the estimate seems to work fairly 
well. 

One remarkable thing is that in our observations a 
granular eruption (almost) never coincides with a thick- 



In conclusion, we have studied the influence of the 
boundaries on the different phenomena that can be ob- 
served after impact of a ball into a loosely packed sand 
bed: The penetration of the ball into the bed, the forma- 
tion of a void, its collapse and the creation of a granular 
jet, the shape of the granular jet, and the presence of 
a granular eruption. We have shown that all of the ob- 
served behavior is generally well captured by the drag 
law and hydrostatic collapse model of Section [Til 

More detailedly, in the first part of this study, we have 
shown that the proximity of the bottom changes these 
phenomena, starting with the obvious modification of the 
final position of the ball, which below a certain depth just 
hits the bottom. The height of the jet is affected, when 
the void closure is constrained to happen on top of the 
ball. A granular eruption at the surface only happens if 
the volume of the entrapped air bubble is large enough, 
and can be fully suppressed by decreasing the height of 
the sand bed. 

In the second part we have investigated the influence 
of nearby side walls. Here we find a strong influence on 
the drag force that the sand bed exerts on the ball when 
it moves through the sand bed: We find that the hy- 
drostatic drag force component becomes less important, 
whereas the quadratic (velocity-dependent) component 
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becomes more important. The latter can be traced back 
to the increased importance of the air flow in the con- 
tainer due to the confinement. Apart from the question 
why and how the coefficients depend on ambient pres- 
sure and container diameter, the drag model of Section HI1 
provides a quite accurate description of the observations 
for most of the parameter space. Only the results for 
the smallest container at low ambient pressure cannot 
be explained using this framework, due to the constant 
velocity plateau that is observed during the motion. 

The formation and subsequent collapse of the cavity 
is not only influenced by the modification of the trajec- 
tory of the ball; also a smaller amount of sand is involved 
in its collapse which therefore takes longer for decreased 
container size. Apart from this, the simple hydrostatic 
collapse model of a cylindrical cavity presented in Sec- 
tion [IT] accounts well for most of the observations. In this 
way, the modification of the closure time, and closure 
depth observed in our experiments, can be understood. 

As a result of both the changes in the ball's trajectory 



and the smaller amount of sand that is involved in the 
collapse, the jet height is affected by the proximity of 
the wall. In the parameter range of our experiments the 
unconfined behavior is retrieved when the diameter of the 
container is larger than 7d; this value however does seem 
to depend on the Froude number, and is larger when the 
Froudc number is larger. The occurrence of a granular 
eruption was shown to be correlated with the size of the 
air bubble entrapped inside the sand bed. 

Finally, this paper culminates in the proposal of a new 
mechanism for the formation of the thick-thin structure, 
based upon a second collapse that occurs on top of the 
ball when it has come to a standstill. To obtain a thick 
part in the jet, the second jet coming from this secondary 
collapse needs to be formed fast enough to penetrate the 
rapidly growing layer of sand that is being created around 
the point where the first jet had originated. 

The work is part of the research program of FOM, 
which is financially supported by NWO; S. v. K. and S. 
J. acknowledge financial support. 



[1] H.M. Jaeger, S.R. Nagel and R.P. Behringer, Rev. mod. 

Phys., 68, 1259-1273 (1996) 
[2] ST. Thoroddsen and A.Q. Shen, Phys. Fluids 13, 4 

(2001) 

[3] D.Lohse, R. Bergmann, R. Mikkelsen, C. Zeilstra, D. van 
der Meer, M. Versluis, K. van der Weele, M. van der Hoef, 
H. Kuipers, Phys. Rev. Lett. 93, 198003 (2004) 

[4] J.R.Royer, E.I. Corwin, A. Flior, M.-L. Cordero, M.L. 
Rivers, P.J. Eng, H.M. Jaeger, Nature Phys. 1, 164 
(2005) 

[5] J.R.Royer, E.I. Corwin, P.J. Eng, H.M. Jaeger, Phys. 

Rev. Lett. 99, 038003 (2007) 
[6] J.R.Royer, E.I. Corwin, B. Conyers, A. Flior, M.L. 

Rivers, P.J. Eng, H.M. Jaeger, Phys. Rev. E 78, 011305 

(2008) 

[7] G.Caballero, R. Bergmann, D. van der Meer, A. Pros- 
peretti, D. Lohse, Phys. Rev. Lett. 99, 018001 (2007) 

[8] J.O. Marston, J.P.K. Seville, Y.-V. Cheun, A. Ingram, 
S.P. Decent, M.J.H. Simmons, Physics of Fluids 20, 
023301 (2008) 

[9] J.S. Uehara, M.A. Ambroso, R.P. Ojha, D.J. Durian, 

Phys. Rev. Lett. 90, 194301 (2003) 
[10] M.P. Ciamarra, A.H. Lara, A.T. Lee, D.I. Goldman, I. 

Vishik, H.L. Swinney, Phys. Rev. Lett. 92, 194301 (2004) 
[11] D.Lohse, D. Rauhe, R.P.H.M. Bergmann, D. van der 

Meer, Nature (London) 432, 689 (2004) 
[12] J. R. de Bruyn and A. Walsh, Can. J. Phys. 82, 439 

(2004) 

[13] M. Hou, Z. Peng, R. Liu, K. Lu, C.K. Chan, Phys. Rev. 

E 72, 062301 (2005) 
[14] L.S. Tsimring and D. Volfson, in Powders and Grains 

2005, edited by R. Garcia-Rojo, H.J. Herrmann, and 

S. McNamara (Taylor and Francis, London, 2005), p. 

12151223 

[15] Hiroaki Katsuragi and Douglas J. Durian, Nature Physics 
3, 420-3 (2007) 

[16] E. L. Nelson, H. Katsuragi, P. Mayor, and D. J. Durian, 
Phys. Rev. Lett. 101, 068001 (2008) 



[17] A. Seguin, Y. Bertho, and P. Gondret, Phys. Rev. E, 78, 

010301 (R) (2008) 
[18] S.J. de Vet and J.R. de Bruyn, Phys. Rev. E, 76, 041306 

(2007) 

[19] R. Mikkelsen, M. Versluis, G.W. Bruggert, E.G.C. 

Koene, D. van der Meer, K. van der Weele, D. Lohse, 

Phys. Fluids, 14, S14 (2002) 
[20] D. Gilbarg and R.A. Anderson, J. Appl. Phys. 19, 127 

(1948) 

[21] R.P.H.M. Bergmann, D. van der Meer, S. Gekle, A. van 
der Bos, and D. Lohse, J. Fluid Mech., in press (2009) 

[22] Note that the quadratic drag is called "inertial drag" and 
a = m/do where m is the mass of the sphere and do is 
the constant introduced in [l5| . 

[23] Starting from the collapse time t c , the rise time of the 
bubble has been estimated as that of similarly sized bub- 
ble in a liquid experiencing Stokes drag, assuming that it 
rises in a straight path with its terminal velocity imme- 
diately, i.e., from the balance 4>PgVbubbi c g oc ^V^bie'W 
(with the packing fraction (j) and the dynamic viscosity 

n assumed to be constant) we have ^risc OC Fj-j^k^^p with 
Fbubbic oc (z(t c ) — z c )d 2 . Now, we estimating the initial 
position of the top of the bubble as zf — (z(t c ) — z c ). 
This leads to t\ — (zs + z c — z(t c )) /v r i sc oc (zf + z c — 
z(t c ))(z(t c ) — z c )~ 2/3 d~ 4/3 g~ 1 . The proportionality con- 
stant was fitted to give the correct large depth behavior. 
Regarding the dissolution time up to hhed ~ lid we 
can estimate the pressure difference by the hydrostatic 
pressure in the center of the bubble at the moment the 
ball has stopped, i.e. AP « (f>p g g(z{ — (l/2)(z(t c ) — z c )) 
(where z(t c ) equals h bc d — d for h bc d < 5.5 d), the volume 
of entrapped air again as Fbubbic oc ( z (t c ) — z c )d 2 , and, 
since ball reaches the bottom, the shortest path is around 
the ball through the sintered plate, i.e., H ~ d. Using 
Darcy's law we have t 2 « Vbubbic/Q oc Vbubbic-ff/AP. 
Inserting the above quantities we obtain ti oc (z(t c ) — 
z c )d 3 /{(j)p g g(zf - (l/2)0(t c ) - z c )). Above h be d = lid 



15 



only the path length changes to H oc hhcd — Zf such that 
fa oc (z(t c )- z c )(h bcd - z f )d 2 / (<j>p g g(z{- (l/2)(z(t c )- z c )). 
Again, the proportionality constant was used as a fitting 
parameter. 

[24] The (small) differences in the fitting parameters n and 
a found for the various Froude numbers were consistent 
with the measurement error, except for the smallest con- 
tainer diameter at the smallest pressure, as explained in 
the text. 

[25] To obtain Fig. [14] the closure depths z c for the diameter 
of D = 8.5 cm (which were not measured directly) are 
obtained by interpolation from Fig. Ilia . 

[26] In Fig. [14] we choose to plot the final depth zi on the 
vertical axis rather than the ambient pressure p to sup- 
port the analysis in this Section which makes use of the 



final depth. Since, as we have shown earlier, zt mono- 
tonically decreases with decreasing pressure, the phase 
diagrams look very similar when using p instead of Zf as 
the quantity plotted on the vertical axis. 

[27] Note that the largest container size (the one without an 
inserted cylinder) has not been included because of its 
square cross section, which is found to have a marked 
influence on the jet shape. 

[28] At this point it is good to note that such a mechanism ex- 
plains why the occurrence of a thick-thin structure never 
seems to interfere with the jet height: The height is de- 
termined by the free flight of the thin part which is being 
formed at the first closure z c , i.e., before the formation 
of the thick part can become of influence. 



